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Introduction 

In the present investigation we aim to provide experimental information on and thereby 
understanding of the generation and propagation of spark-ignited, outwardly propagating cellular 
flames, with three major focuses. The first is to unambiguously demonstrate the influence of the 
four most important parameters in inducing hydrodynamic and diffusional-thcrmal cellularities, 
namely thermal expansion, flame thickness, non-unity Lewis number, and global activation 
energy. The second is to investigate the critical state for the onset of cellularity for the stretch- 
affected. expanding flame (1-41. The third is to identify and consequently quantify the 
phenomena of self-acceleration and possibly auto-turbulization of cellular flames |5. 6|. Due to 
space limitation the effects of activation energy and the critical state for the onset of cellularity 
will not be discussed herein. Details can be found in Ref. 7. 

Experiments were conducted using CjHg-air and fL-Oj-N’ mixtures for their opposite 
influences of nonequidiffusivity. The additional system parameters varied were the chamber 
pressure ip) and the mixture composition including the equivalence ratio (f>). From a sequence of 
the flame images we can assess the propensity of cell formation, and determine the instantaneous 
flame radius <R). the flame propagation rate, the global stretch rate experienced by the flame, the 
critical flame radius at which cells start to grow, and the average cell size. 

Effects of Instability Parameters 

The hydrodynamic theory of Darricus |8| and Landau |9) shows that, in the limit of an infinitely 
thin flame propagating with a constant velocity, the flame is unstable to disturbances of all 
wavelengths. The growth rate is proportional to the density jump across the flame, increasing 
with increasing <T. Thus a is probably the most sensitive parameter controlling the onset of 
hydrodynamic instability. 

Next to o. the flame thickness S_ \ is also expected to have a strong influence on the 
hydrodynamic instability, for two reasons. First, it measures the influence of curvature which, 
being positive for the expanding flame, has a stabilizing effect on the cellular development. The 
thinner the flame, the weaker is the influence of curvature and consequently the stronger is the 
destabilizing propensity. The second influence is that it controls the intensity of the baroclinic 
torque developed over a slightly wrinkled flame surface, which depends on the density gradient 
across the flame and the pressure gradient along the flame (10J. Since the density gradient 
increases with decreasing flame thickness, development of the hydrodynamic instability is 
correspondingly enhanced due to the increased intensity of the induced baroclinic torque. 

For the development of the diffusional-thcrmal instability, an appropriate parameter 
representing the effect of nonequidiffusion is the flame Lewis number. Le [11|. It is well 
established and understood theoretically that unstrctchcd flames arc diffusionally unstable 
(stable) for Le ' s that arc smaller (greater) than a value slightly less than unity. 
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The last parameter of importance, especially for the present outwardly expanding flame, 
is the Karlovitz number, which is the nondimensional stretch rate and is defined as Ka = 
(2/R)(dR/dt)/(s u "/Sr). where s u ° is unstrctchcd laminar burning velocity and / is time. It was 
theoretically shown for the stagnation flame |I2) and the expanding flame |I. 2] that the 
associated positive stretch tends to be stabilizing. Conceptually, cells cannot form if their growth 
rate is smaller than that of flame expansion. Since the expanding flame suffers the strongest 
stretch during the initial phase of its propagation when its radius is small, the tendency for cell 
development is expected to increase as the flame propagates outward. The influence of stretch 
also shifts the critical Le. at which the flame response reverses, to a smaller value than that of the 
unstrctchcd flame. 

We now present in sequence the influence of Sj. a. and Le on the propensity of cell 
development, with the comparison conducted at the same or similar Karlovitz numbers. 

FUunc Thickness Effects : Figure 1 shows the burning sequences of almost diffusionally neutral, 
stoichiometric propane-air flames in 2. 5. and 10 atm. The corresponding flame thicknesses arc 
Sr = 0.17. 0.10. and 0.062 mm. The sequence at 2 atm. shows that the flame surface remains 
smooth, with the presence of only a large ridge throughout the observation period. At 5 atm. a 
few large cracks, formed as a consequence of the disturbance caused by the spark discharge, 
persist down to Ka = 0.05. However, when Ka is reduced to 0.03. small cells of an average size 
of 3.6 mm emerge. The sequence at 10 atm. shows that the cracks first grow in a self-similar 
manner up to Ka = 0.08. Then further cracking through branching is developed, as seen for Ka = 
0.05. Eventually the average cell size is reduced to 1.1 mm at Ka = 0.03. The above flame 
morphology readily substantiates the concept that positive stretch and a thicker flame tend to 
delay the onset and development of hydrodynamic cells, and that as the flame becomes thinner 
not only it is destabilized earlier (i.c. for larger values of Ka). but the cell size is also smaller. 

Thermal Expansion Effects: Figure 2(a) shows the burning sequences of two diffusionally stable 
flames that have different values of <rbut nearly the same values for other parameters: namely a 
0 = 0.9. p = 10 atm., propane-air flame with <7 = 7.7. S t = 0.068 mm. and Le = 1.6. and a <p= 3.0. 
p = 5 atm., hydrogen-air flame with <7= 5.9. Sr = 0.059 mm. and Le = 1.7. Since the former has a 
larger expansion ratio, it is more unstable, as shown. 

Nonequidiffusive Effects: For this demonstration, we compare in Fig. 2(b) the burning sequences 
of two p = 5 atm.. <7= 7.7. propane-air flames at 0= 0.9 (Le =1.6 .Sr= 0.10 mm) and 1.5 (Le = 
0.95. Sr = 0.21 mm), such that the former is diffusionally stable while the latter unstable. 
Considering that the former has a smaller Si and hence is hydrodynamically more unstable, the 
fact that the latter actually exhibits a more prominent instability pattern illustrates the powerful 
destabilizing influence imposed by the diffusional-thermal instability mechanism. The absence 
of any cells for the diffusionally stable flame, even at Ka = 0.10. and the presence of the (large) 
hydrodynamic cells for the diffusionally unstable flame at Ka = 0.32. possibly triggered by the 
diffusional-thermal instability, arc particularly worth noting. 

Self-Acceleration of Cellular Flames 

A series of crucial questions can be asked about the propagation rate of the expanding flame. 
First, after development of the cells, will the wrinkled flame propagate faster than the original 
smooth flame? Second, if it is faster, will it accelerate? Third, if it accelerates, will the 
acceleration be constant? Stated alternatively, by determining the flame radius history R(t) - t". 
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and hence the propagation velocity dR(t)/dt ~ affirmative answers to the above three 
questions would require a assuming a constant value greater than unity. Furthermore, a constant 
a would imply that the flame propagation and morphology have a fractal character, and that if 
the fractal dimension is close to that of turbulent flame propagation, the self-acceleration process 
of the wrinkled flame propagation can be considered as one of auto-turbulization |5. 6. 13]. 

Figure 3 plots R(t) for three highly unstable flames of H2-15%0’-N2 mixtures of 0 = 0.5. 
1.0. and 1.5, at 15 atm. The corresponding Le arc 0.40. 0.94. and 1.65. The data were 
subsequently fitted according to R(l) = R a + At", where R„ is the virtual origin, representing the 
effect of the initial period of steady propagation. The fitting yields R = 0.0013+1. 90/' 26 for 0 = 

0. 5. R = -0.0015+41.53/' 36 for 0= 1.0: and R = -0.0032+26.04/ 1 2 ‘ for 0= 1.5. where R is in m 
and / in see. The correlation of the data in the numerical fitting is accurate to 99.9%. The inset 
shows the corresponding Rlt)-R„ in logarithm coordinates, demonstrating the linearity of the 
fitting and hence the constancy of a. 

The above fitting shows that all the three flames arc self-accelerating. The acceleration 
exponents determined for the three eases arc 1.26. 1.36. and 1.23, which do not show any 
particular pattern in relation to any of the flame parameters. Indeed, the three values arc 
purposely shown as typical values since at the present state it is difficult to pin down a precise 
number, if indeed there exists a precise number. However, they are smaller than the value of 1 .5 
reported in Ref. |5|. but arc basically within the range of 1.25 to 1.5 mentioned in a subsequent 
publication 1 14]. Thus by using the relation a= 1/(1 -d) |6], where d is the fractal excess, our 
results basically yield d in the range of 0.20 to 0.25. implying a fractal dimension. (2 +d). of 2.20 
to 2.25. 

Summary 

From high-speed imaging of outwardly propagating propane-air and hvdrogcn-oxygcn-nitrogcn 
flames under elevated pressures, it was demonstrated that hydrodynamic instability is greatly 
enhanced with increasing pressure and hence decreasing flame thickness. The cellular flames 
were found to be self-accelerating, including those that arc diffusionally unstable, with fractal 
dimensions between 2.20 and 2.25. 
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Figure I . Schlieren photographs of stoichio- Figure 2. Schlicrcn photographs of (a) CjHg- 
mctric C.iHs-air Haines at 2. 5. and 10 atm. air flames of 6 = 0.9 at 10 atm. and H 2 -air 














